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We study the UV-scaling of the flavorless gluon propagator in the Landau gauge in an energy window up to 9 
GeV. Dominant hypercubic lattice artifacts are eliminated. A large set of renormalization schemes is used to test 
asymptotic scaling. We compare with our results obtained directly from the triple gluon vertex. We end-up with 
— 318(12)(5) MeV and 292(5)(15) MeV respectively for these two methods, compatible which each other 
but significantly above the Schrodinger method estimate. 
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1. Definition of the coupling constant from 
the propagator 

The Euclidean two point Green function in mo- 
mentum space writes in the Landau gauge: 

G'^^^'iP, -P) - G(-'\p')Sa, (s,^ - ^) (1) 

where a, b are the color indices. 

In any regularization scheme (lattice, dimen- 
sional regularization, etc.) with a cut-off A (a~^, 
(d — 4)~^, ...) the bare gluon propagator in the 
Landau gauge is such that 



r(p2) = lim 



(2) 



is independent of the regularization scheme. 

Lattice calculations provide us with a measure- 
ment of the bare propagator, and hence of the 
observable T in eq. (^). 

We define the running coupling constant, 
a(p^), through: 
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71 



(3) 



An (47r)2 (47r)3 
for a given set of real coefficients 70 , 71 and 72 . 
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Note that in the MOM scheme, the gluon renor- 
malization constant Z3 is given by ^3(9^, A) = 
q^Gbare{q^,^)- Our observable F is then the 
anomalous dimension of the gluon renormaliza- 
tion constant in the MOM scheme. But it is im- 
portant to realize that this does not imply that 
we are stuck to the MOM scheme : F is defined in 
a scheme independent way by eq.(^ and can be 
studied in any general scheme defined by eq.(||). 
Nevertheless this imposes the value of 70 to be 
equal to 13/2 , the first (universal) term of the 
anomalous dimension. 

Our strategy to study 5(/i) on the lattice from 
the propagator is the following: 
• Measure on the lattice the bare gluon prop- 
agator in Landau gauge, extract Gbarov 
following eq.(|l|) and compute ln[Z3(p, a)] 



GlZ^ip,a) 
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• Choose 71, 72 and integrate the two coupled 
differential equations: 

dlnZ3{q,a) 



dlnq^ 



70 ~ 
= a 
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71 
(4^) 



72 



5^(4) 
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-0(55)(5) 



(/32 is known for any 71 and 72 [|l[ ) • The solu- 
tion for Z^{q,a) and a{q) depends on the val- 
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ues Z^{^, a) and a(/i) given at some initial value 

• Fit the lattice data with the previous solutions 
to determine the "best" Z3{fj,,a) and a{^). Ex- 
tract A (we use three-loop expressions here) and 

71 -~i 



270 /3o 



convert to Ajjg through = A exp 

where 7^ ~ 155.3 is the value of the coefficient 71 
in the MS scheme. 



2. Hypercubic artifacts 

In an hypercubic volume the momenta are the 
discrete sets of vectors — where the 

are integer. For a given orbit of the continuum 
isometry group (characterized by a given value 
for n^), there are distinct orbits of the hypercubic 
isometry group. For example — (1,1,1,1) and 
= (2,0,0,0). 

A usual approach to reduce these hypercu- 
bic effects is to use only the "democratic" sets 
of momenta, i.e. those for which the momen- 
tum is equally distributed on the different di- 
rections and/or the use of the lattice momenta 
Pfj^ = ^ sin (-f^). We propose another approach, 
thanks to our large statistics, which allows us to 
eliminate these hypercubic artifacts. On the lat- 
tice, an invariant scalar form factor like G^^^ is in- 



deed a function of the 4 invariants : 



pi' 



with m = 2, 4, 6, 8. We will neglect the invariants 
with degree higher than 4 since they vanish at 
least as a^. 

When several orbits exist for one , we have 
found a nice linear behavior cj^^ (p^, a^pW) = 

It is then 



GL1(^'^o) 



dG 



possible to extrapolate to p^'*' — and we can 



ift (P^ ' -P'''' ) iu all our 
analysis. We have introduced this method in H]. 
To reach higher energies and exploit also the cases 
where only one orbit is available, the method has 
been extended : from dimensional argument, the 
slope ^23p[4 i is expected to behave like 1/p^ as 



L 

ansatz 



g, we found that the 
^(1 + c exp(-dLp)) 



CXI. As shown in 

OG ^ 

gives a good description at fixed (3 of all our data 
thus allowing us to extrapolate to p^^' — 0. 



3. Finite volume effects 

In 1^, with our lattices a,t (i — 6.2 and 6.4, 
we found that the gluon propagator has not 
yet reach asymptotic scaUng at three-loop 
even at yu ~ 5 GeV. Higher energies are then 
required to test scaling. We have now lattices at 
(3 = 6.8, but as we are limited to 24^ which is not 
a too large volume at this value of /3, attention 
has to be paid to the finite volume effects. 
To that end we have generated lattices at (3 — 6.0 
with several volumes 12^, 16^^, 24''^ and 32'* and 
tried several functional forms to parametrize the 
finite volume effects at fixed value of (3. Finally 
we proposed: 



G{ft^(p',oo,a) 



1 



(6) 



V2 e 



-D3 Lp 



which gives a good parametrization of all the (3 = 
6.0 data in the UV part. See figure (|l]). 

Once the parameters vi,V2, W3 have been de- 
termined from the lattices at /? = 6.0, the 
parametrization (^) can be applied to extrapo- 
late our results on 16** and 24** lattices at p = 6.8. 
The good agreement shown by the curves result- 
ing from the extrapolation (figures can be found 
in ^) is a crucial test for the validity of such a 
parametrization for the finite volume effects. 

4. Results 

To calibrate the lattice a.t f3 — 6.2, we have used 
the value for the the lattice spacing which has 
been measured recently with a non-perturbatively 
improved action (free from 0(a) artifacts) 
The measured value is : a-i(/3 = 6.2) = 2.72(11) 
GeV. Other lattices have been calibrated rela- 
tively to the one at /3 = 6.2 with the results 
for a^/a, the string tension in lattice units, pub- 
lished in §. We took: a-^{(3 = 6.0) 1.94 GeV, 
a'^{l3 = 6.2) = 2.72 GeV , a-i(/3 = 6.4) = 3.62 
GeV and a-\P = 6.8) = 6.03 GeV. 

We have applied the strategy described above 
to study a, using the data at P = 6.8 and found 
at three-loop using the MOM scheme: a}-^"^'' ~ 

MS 



0.346 A^~^^ = 315 ± 12 MeV, 

MOM 

only statistical here. 



where the error is 
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We can exploit the large scheme space to which 
we have easily access to define a domain of "good 
schemes" by imposing some constraints on 71 and 
72- The idea is to select schemes for which the 
successive terms in the expansions in a for the 
/3-function and F are not exceedingly large com- 
pared to the previous terms (see [|| for the precise 
definition of the cuts we used). A large scheme 
dependence in this set of schemes would indicate 
that we are still far from perturbative scaling. We 
found that mom is a "good scheme" while ms is 
not and the result for j^^^°p^ ranges from 313 
to 323 MeV in the domain oP'good schemes" . We 
consider this dispersion to give an estimate of the 
effect of higher order terms in the perturbative 
expansions and the fiavorless extracted from 
the propagator up to 9 GeV we quote is then: 

-^'^^■'^SIS (12) (5) 



-MeV (7) 



6.03 GeV 

We can compare this result to the one we obtain 
by direct standard textbook methods from the 
triple gluon vertex in the mom scheme j^]. Data 
for a are less precise than those for the propa- 
gator so we have not yet tried to eliminate the 
finite volume effects with a method similar to the 
one used here. Consequently we did not use our 
/3 = 6.8 lattice and stay with our large volume 
lattices at smaller (3. A was found to be nearly 
independent of /? for (3 ~ 6.0, 6.2 and 6.4 |jl|] : 

a-i(/3 = 6 



A 



(vertex) 



292 (5) (15) MeV- 



(8) 



MS V / V / g_Q3 Q^y 

Our two methods, based on Green functions, give 
compatible results, both significantly larger than 
the one obtained over a large energy range with 
the Schrodinger functional ||6|]. 
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Figure 1. Plot (a) shows the free-hypercubic prop- 
agators evaluated on 12^ (black squares), 16'* 
(black circles), 24* (white circles), 32* (white 
squares) lattices at (3 — 6.0. Plot (b) shows the ex- 
trapolation L — !■ 00 0/ these free-hypercubic propa- 
gators using the parametrization given in Eq. (j^. 



